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Abstract 

We consider a moving interface that is coupled to an elliptic equation 
in a heterogeneous medium. The problem is motivated by the study of 

00 displacive solid-solid phase transformations. We show that a nearly flat in- 

terface is given by the graph of the function g which evolves according to the 
equation gt(x) = —{—/S) 1 ^ 2 g{x) + ip{x,g(x)) + F. This equation also arises 
in the study of dislocations and fracture. We show in the periodic setting 
that such interfaces exhibit a stick-slip behavior associated with pinning 

CN and depinning. Further, we present some numerical evidence that the effec- 

tive velocity of the phase boundary scales as the square-root of the excess 
macroscopic force above the depinning transition. 



1 Introduction 



Hysteresis is ubiquitous in materials science, and is associated with nucleation 
and propagation of interfaces and defects. This paper concerns the propagation of 
interfaces immersed in an elastic medium. In particular, the paper is motivated 
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Figure 1: Pinning-depinning behavior as calculated from a one-dimensional model 
in [2]. The interface is stuck up to a critical force F*, and breaks free above it 
with a particular scaling. For large F, the average velocity is linear in the applied 
force. 

by phase boundaries in solids that undergo a displacive phase transformation such 
as the martensitic phase transformation. In these transformations, one has phase 
boundaries across which the crystal structure changes without any diffusion or loss 
of compatibility. Many interesting properties of such materials, like the shape- 
memory effect, are associated with the nucleation and evolution of these phase 
boundaries. As the interface propagates, the change in crystal structure potentially 
gives rise to elastic fields. Thus, one has a moving interface problem that is 
coupled to an elasticity problem. Further, every material contains defects like non- 
transforming precipitates which makes the medium inhomogeneous. An important 
question is the role of these defects, and this motives our current work. Similar 
issues arise in ferroelectricity, ferromagnetism, and other phenomena. 

There is a well-developed framework to study the evolution of martensitic phase 
boundaries, and this is described in detail in the recent monograph of Abeyaratne 
and Knowles [3 J . Briefly, one defines a thermodynamic driving force either through 
the rate of dissipation or through the variation in the total energy with respect 
to the position of the interface, and then postulates a kinetic relation that relates 
the driving force to the normal velocity of the interface. Microscopic theories 
suggest that the kinetic relation has viscous character passing smoothly through 
the origin [21 121] • Such a kinetic relation predicts that the hysteresis goes to zero as 
the rate of loading goes to zero. However, experiments clearly show otherwise: the 
hysteresis does not go to zero with loading rate and instead settles on a non-zero 
value independent of loading-rate for slow enough rates. Such observations suggest 
a stick-slip behavior where the interface is stationary below a critical driving force 
and moves freely above it. It is often suggested that pinning of the phase boundary 
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by defects is responsible for this transition from microscopic viscous to macroscopic 
stick-slip behavior. 

A one- dimensional calculation, as found in [U [5 J illustrates how a local wiggly 
potential can pin a phase boundary and lead from a linear kinetic relation to 
a stick-slip behavior. Assume a bar with a 1-periodic local driving force <p(x) 
(smooth and with non-degenerate global maximum and minimum), and assume 
that the velocity of the interface is given as v = ip + F, where F is the constant 
external applied force. The amount of time it takes for the interface to travel one 
period can now easily be calculated to be 



if F > — mimp or F < - max ^. Otherwise (i.e., if — mm tp < F < — max tp), the 
time is infinite and the interface is stuck. Further, close to the critical F, say F pa 
— max (p, the interface is slow only in a few isolated points but propagating freely 
everywhere. This implies, under some non- degeneracy and regularity conditions, 
that the effective velocity scales as the square-root of the excess force. Thus, 
the effective velocity v — ^ of the interface now exhibits a behavior of the form 
shown in Figure [T] A rigorous proof of the transition from a viscous microscopic 
kinetic law to a rate independent evolution through the interaction with a wiggly 
potential can be found in [19] . The question whether such a stick-slip behavior is 
also observed in models for phase transformations in higher dimension motivates 
this work. 

We present a sharp interface model for the quasistatic evolution of a martensitic 
phase boundary in higher dimensions in Section 2. We limit ourselves to the scalar 
anti-plane shear setting (where the displacements are scalars) though the ideas 
and results hold for the general case. In this model, a free boundary separates 
two material phases. Each phase is characterized by a distinct transformation 
or stress-free strain where the elastic energy density admits its minimum. We 
also assume that the material contains a number of non-transforming precipitates. 
Importantly, both the phases as well as the non-transforming precipitates have the 
same energy. A similar model was studied by Craciun [6] . 

We then derive an approximate model for a nearly flat interface. We show, 
using methods of T-convergence, that the elastic energy of a nearly flat interface is 
approximated by the if 1//2 -norm of a function whose graph describes the interface 
(Theorem 2.1). We also argue that at low volume fraction, the precipitates give 
rise to a local forcing which scales similarly to the elastic energy. We thus conclude 
that the interface is described by the graph of a function g which is governed by 
the equation 





g t = -(-A) 1 / 2 g + i P (x 1 



...,x n ,g(xt, ...,£„)) + F 



(2) 
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for a given (p : W i+1 — > M with zero mean. On the periodic domain we consider, 
this equation may be compactly written by its Fourier series, 

g t (k) = -\k\g(k) + 0(k)+F{k). (3) 

From now on, g indicates the Fourier series of the periodic function g. Note that 
the equation is still nonlinear, since the driving force ip depends on g. 

While we derive this model from phase transformations, it has also been used 
to study dislocations [18] (see [H] for some rigorous analysis of that model) as well 
as fracture [20]. In fact, a very similar model is derived in [13J as the homogenized 
limit of an interacting system of individual dislocations. Front-type solutions as 
limits of reaction-fractional-diffusion equations are derived in [T6] . A closely related 
parabolic model, 

g t = Ag + <p + F (4) 

has been used to study pinning of surface energy dominant interfaces by defects. 
The large physics literature has concentrated on the situation where (p is random, 
and has shown using scaling arguments and numerical simulation that these equa- 
tions lead to a pinning/ depinning transition with a critical exponent which varies 
from situation to situation [1]. Dirr and Yip |9] presented a rigorous analysis of 
the parabolic model (El) in the periodic setting (99 is periodic). A rigorous analysis 
of the random case remains the topic of ongoing research (see for example pj [S] ) . 

We study the behavior of the solutions of ^ in Section 3 where cp is (1— )periodic 
following Dirr and Yip [9]. We show that there is a critical F* > such that ^ ad- 
mits a stationary solution for all F < F* (Theorem 3.5). Further, for each F > F*, 
there exists an unique T such that ^ admits a space-time periodic solution (The- 
orem 3.7 and Proposition 3.9). Thus, we may regard 1/T as the effective velocity 
of the interface. In Section 4, we discuss the behavior of the effective velocity near 
the depinning transition and present some numerical examples indicating that the 
effective velocity scales as the square-root of the excess force F — F* in the case 
of a smooth, non-degenerate heterogeneity. 

2 A model of phase transformations 

2.1 Phase transformations in the presence of defects 

We consider a model proposed by Craciun [6]. Since we are interested in the over- 
all propagation, we consider the domain to be a strip, O = T" x M where T n is an 
n-dimensional torus as shown in Figure [2] The domain Q is divided into two parts, 
E and Q\E occupied by two phases, and separated by the phase boundary T of 
codimension 1. The domain also contains a number of non-transforming precipi- 
tates, occupying the set Aj. The two phases are characterized by two stress- free 
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Figure 2: A phase boundary in a strip containing non- transforming precipitates. 



strains £ e M. n+ , and the non-transforming precipitates are characterized by the 
stress-free strain £°. We further assume that all phases and precipitates have equal 
elastic modulus (which we take to be identity without loss of generality). Thus 
the elastic energy of domain subjected to the displacement u : Q — > R is given by 

Elastic = / \ 1 1 Vw - £ E \ 1 2 &xdy. (5) 

where 

( £+ (x,y)eE\\J i A i 

£e{x,v) = { r (x,y)e(Q\E)\[j i A i (6) 

{ e (x,y)e\J i A i . 

and V denotes gradient with respect to (x,y). 

For a given interface T, we obtain the displacement by minimizing the energy 
subject to appropriate boundary conditions. This displacement satisfies the Euler- 
Lagrange equation 

Au = div£ E . (7) 

We then say that the interface T evolves in a specified manner that reduces the 
(optimal) elastic energy [3]. 

In this section, we seek to find an approximation for this particular model. To 
motivate the approximation, consider the situation where one does not have any 
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non-transforming precipitates and the interface that minimizes the total energy. 
It is easy to show in the current scalar setting that the interface is a plane with 
normal n = (£ + — — £~| Q Therefore, we expect that an arbitrary interface 

will soon become almost planar with this normal and then evolve in an almost 
planar manner. Further, in the presence of defects, we expect the interface to be 
distorted close to them due to the elastic fields created by the defects but be largely 
planar away from them. Thus, if the concentration of defects is small, the interface 
remains largely planar. This is consistent with numerical observations [10J. All 
of this motivates us to seek an approximation for the model in the case that the 
interface is almost planar. 

Due to the linearity of the Euler-Lagrange equation ([7]), one can split the 
transformation strain and the displacement into components depending only on 
the interface and on the precipitates, respectively. Take 



u = u p + u p ; ^ E = e + e 



and fix 



and 



~ { |- xe(Q\E), (9) 

The displacements u p , u r solve the Euler-Lagrange equations associated with £ r , £ p , 

Au r = div£ r ; Au p = div£ p . (11) 

This fixes the functions u v and u p up to an affine component. 

Substituting this decomposition back into the elastic energy, we can expand it 
as follows: 

'"IVu-^H 2 = / i||V M r -e r || 2 (12) 



n 2 2 



+ / (W-O-v?/ 

Jn 

- [ - e) ■ ?■ 
h\vu p -e\\ 2 



The first term in ( 12 ) is the energy associated with the interface T in the absence of 



any precipitates. We shall call this the self-energy of the interface, and we find an 



"hi the vectorial setting of linearized elasticity, there are two possible normals [17] 
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approximation for it in section 2.2 We can use the Euler-Lagrange equation ( 1 1 ) to 
show that the second term is identically zero, assuming that the affine component 
of u v is chosen appropriately. The remaining terms describe the interaction of 
the precipitates with the interface, which is dealt with in a heuristic manner in 
section 12.31 

We seek approximations to the self energy and the interaction energy when the 
interface is almost planar. We specialize to the case when 

e ± = (0,...,0,±l/2) T (13) 

and the preferred normal is n = (0, . . . , 0, 1) T . We assume that the phase boundary 
is the graph of a H l and essentially bounded function g, i.e., 

T = {(21,22, .-,x n ,y) such that y = g{x u ...,2 n )} , (14) 

and the transformed domain E = {(x,y) : y < g(x)}. 



2.2 Approximation of the self-energy 



We first show that the self-energy, the first term on the right hand side of (12), 
scales as e 2 if the interface is rescaled by a factor e and in the limit may be ap- 
proximated by the one half the square of the if 1 / 2 semi-norm of g. Since this term 
is independent of the precipitates, we ignore these for now. It is convenient to 
subtract the piecewise constant, curl-free, function |(1 — ff(?/))jT]from the transfor- 
mation strain in order to make it supported in a bounded region around y — 0. 
This does not change the energy of minimizers of the energy functional, since one 
can simply subtract the respective integral from u. From now on, we thus assume 
that, for a given function g: T n — > R, 



( \ 



for y < and y < g 
— 1 for y < and y > g 

1 for y > and y < g 
for y > and y > g 



(15) 



and ig{x,y) 



In the following, we take g = eg for g G H (T n ) n 



L°°(T n ), e small, and prove the convergence of the funtional 



1 

^2 



(16) 



*H(y) denotes the Heavyside step function. 
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in the sense of T-convergence with respect to the usual L\ oc metric, to the limit 
functional 

T(u) = l ^\{y=o} '2 1 1 Vm I I' if M = 9 a-e. 
[ oo otherwise. 

Both functionals are assumed to be infinite if u ^ H^ oc . The expression [u] denotes 
the jump of u across {y = 0} (i.e., the difference of the respective traces). We 
write 



ut(x,y) = u £ (x,y + e) 
u £ (x,y) = u £ (x,-y-e) 



Without loss of generality, we assume here that 1 \g\ \ Loo < 1, otherwise the definition 
of uf(x,y) needs to be adapted so that the interface does not penetrate outside 
the cutout region. 

Theorem 2.1 (Approximation of the energy). Let (ej)j e ^ be a decreasing sequence 
of positive real numbers converging to zero. The following assertions hold. 

i) Consider functions u £ . such that J r £j {u £ .) is uniformly bounded. Then there 
exists a subsequence (which we relabel and also index by j) and functions e 
H\ oc (T n x (0,oo)) , such that 

ut — 1 u + 
— 1 u~ 

weakly in Hl oc {T n x (0, oo)) (and thus strongly in L\ oc ), modulo a constant function. 

ii) Given a sequence u £j , such that uf. — > in L\ oc {T n x (0, oo)), we have 

liminf T ei {u £j ) >F(u), (18) 

where u(x,y) = ■u ± (x, ±y) for y > (+) or y < (— ), respectively. 

Hi) Given u e iJ 1 (SI \ {y — 0}), there exists a sequence u £ . so that uf. — >■ u in 

L\ oc {T n x (0, oo)), u- ^ u(x, -y) m L\ oc {T n x (0, oo)), and 

limsup F ej {u £j ) < F{u). (19) 

Proof, i) Clearly, J T n x ( 0oo ) ||Vm^|| 2 is uniformly bounded. This yields the desired 
compactness. 

ii) If [m] = g, the result is immediate from the definition of the energies. What 
remains to show, is that the elastic energy necessarily blows up if the difference of 
the respective traces of u £ (x, e) and u £ (x, —e) does not converge to the function g. 
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Note that, for all e, one can calculate 



> 



£ 1 
£ 1 



> 



> 



T n J-e 
1 

2£ 



VW £ ^ £ g 

e 



d 1 1 

— -u £ - -H(eg(x) -y) + -H(-y) 
ay e e 



1 

2s 
1 

2£ 



^11 

1 

2 



£ 1 
,2 
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— u e - -H(eg(x) -y) + -H(-y) 
ay e e 

d 1 1 
— m £ - -H(eg(x) -y) + -H(-y) 
ay e e 



(20) 



(u e (-,e) -it e (-,-e)) - g\ 



(21) 



Jensen's inequality was used in (20). In (21), u £ (-,y) denotes the trace of u £ on 
(',y). The equality holds, because u e is in H 1 , and thus admits a representative 
that is absolutely continuous on a.e. line [12J. This energy cannot be bounded, 
unless u e (-,e a ) — u £ (-, —e a ) converges to g in L 1 , and thus the jump of u across 
{y = 0} equals g, since weak convergence in H l implies weak convergence of the 
trace. 

Hi) We first assume that the trace of u at y = from below, denoted by T~u, is 
in H 1 . Take 



{x,y) 



u{x,y + e) 

T-u(x) + ^H(eg(x) 

u(x,y-e) 



for y < —e, 
y) for -s < y < e, 
for y > £. 



(22) 



Note that this function is in H^ oc (Q), since the traces at ±e match, inside the 
strip the ^-derivative of u e is bounded by -, and the ^-derivatives of u £ are bounded 
by the sum of the x-gradient of T~u and the x-gradient of g, both of which were 
assumed to be in L 2 . 

The elastic energy J- e (u e ) outside the strip of thickness e remains exactly equal 
to JF{u). The ^/-derivative of the function u £ equals the n-component of £ eg , and 
the x-derivative remains bounded by the (absolute) sum of that of T~u and that 
of g. Thus the integral over the vanishing domain T n x (— e, e) goes to zero. 

In order to obtain the result for arbitrary u G H^ oc , one can employ the usual 
density argument. Functions with iJ 1 -trace are dense and one can approximate 
with the energy bounded by F{u). □ 

The above theorem shows that for a nearly flat interface, the energy due to the 
shape of the phase boundary itself is equal to the energy of a function with a jump 
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of the appropriate height g. It is well known [JS], that the minimum attained at 
u of this energy is equal to one half the H 1 / 2 seminorm squared of g, or, 

min F{u) = \\g\ 2 m/ *. (23) 

My=0=9 

Remark 2.2. Note that we have not proved the T-convergence result for the case 
that g is only in H 1 ^ 2 , since under this weaker assumption a recovery sequence 
can not be found. However, as we will see in Section [3j solutions of the evolution 
problem considered will have the required regularity. 



2.3 Approximation of the interaction energy 

We now turn to the interaction energy, the last two terms of the right hand side 
of @, 

j**:= [ -(yu p -e)-f + h\Vu p -e\\ 2 . (24) 
Jn 2 

Scaling again the transformation strains £ p and £ r by a factor of - and as- 
suming an interface height and precipitate radii rescaled by e, one can see that 
an assumption of a precipitate density of scale e yields an order one term for the 
interaction energy as well. 

Inherently, like the self-energy, J-™* is of course non-local. However, the varia- 
tion of the interaction energy term, as long as the interface does not intersect any 
of the precipitates, is local. Assuming again that the interface is the graph of a 
function g, and £ given as in (13) one can write it as 

6 g r nt = [Vu p (x, g(x)) - f(x, g(x))} • (0, . . . , 0, 1) T =: f(x, g(x)), (25) 

see [3]. From now on we assume that the interaction force for small scarcely 
scattered precipitates can be written in the above local form, ignoring the non- 
locality when passing through inclusions. 



2.4 The approximate model 

The previous sections show that the energy of an interface described by the graph 
of a function g may be approximated as 

£ = \[9\ 2 H ^ + / f" } f(x,y)dxdy. (26) 

^ JT n JO 

The evolution of the interface is described as an L 2 gradient flow of g with 
respect to the energy. Of course, T-convergence of the energy does not neces- 
sarily imply convergence of solutions to the gradient flow, however, we take the 
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T-convergence result in Theorem 2.1 as a good indicator that the evolution of the 
interface can be approximated by a gradient flow of the limit energy. Since the 
variation of the H lj/2 seminorm yields the square root of the Laplacian, we obtain 



g t {x, t) = -(-A) 1/2 g(x, t) + f{x,g{x, t)). (27) 
This equation may be compactly written by its Fourier series, 

g t (k,t) = - \k\g(k,t) + f(^t))(k). (28) 

Note that the equation is still nonlinear, since the driving force / depends on g. 

Finally, we assume that the distribution of precipitates is periodic in the direc- 
tion of propagation of the interface. Thus, the forcing 

f(x, g{x)) = [Vu p (x, g{x)) - £ p (x, g(x))} ■ (0, . . . , 0, 1) T (29) 

can be split up into a periodic term <p(x,g(x)) with zero average and a constant 
term F stemming from boundary conditions at ±oo, so that 

f(x,g(x)) = <p(x ig {x))+F. (30) 



3 Stick-slip behavior 

In this section, we study the stationary equation 

= -(-A) 1 / 2 g + <p + F (31) 

and the evolution equation 

g t = -(-Ay/ 2 g + p + F, (32) 



together with an initial condition, and denote them by {31) p and (32)^, respec- 
tively, indicating the dependence on the behavior on the external force F > 0. 
Our strategy and results closely mirror those of Dirr and Yip [9] who considered 
the Laplacian case. In this section we will make the following assumption on the 
wiggly force ip : T n x R — > K. 

Hypothesis 3.1. The interaction force tp(x,y) is periodic in y, i.e., ip(x,y) = 
<f(x,y + 1), has vanishing average, i.e., J T n x [ 01 ) ^{x.y) = and is Lipschitz con- 
tinuous in both x and y. 
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3.1 Solutions of the evolution equation 

In this section we collect some properties of solutions of the evolution prob- 



lem (32)*-. 



Proposition 3.2. For any initial condition g(-,0) = go G C° , the quasistatic 



evolution problem (32) f admits a unique global in time classical solution. Classical 
solutions to (32 )p admit a comparison principle. Furthermore, for go G L? , a mild 



solution to the evolution problem exists and is classical for any t > 0. 



Proof. We first note that the evolution equation (32) p can be cast in the form 
of [IT] , equation (2), by periodically extending ip to lR n . Since in our model the 
function ip is assumed to be Lipschitz, it follows from [TTJ, Theorem 5, that the 
problem admits a unique Lipschitz continuous viscosity solution, and thus admits 
a comparison principle. We further remark that the fractional Laplacian generates 
an analytic semigroup on the space of continuous functions on T n (for a reference 
see for example [22J, Chapter IX. 11, and note that this property is well known 
for the 'regular' Laplacian), and thus the problem also admits a classical solution 
that can be found via a variation of constants-formula. This also shows that the 
viscosity solution for the problem is periodic. 

For an initial condition only in L 2 , note that the fractional Laplacian also gen- 
erates a semigroup on L p , p > 2, with solutions in H l,p , the domain of —(—A 1 / 2 ) 
for values in LP . We thus can recursively find spaces of higher and higher integra- 
bility for our mild solution at positive time until we obtain a Sobolev-embedding 
into the space of continuous functions. □ 

The next proposition concerns an energy estimate of the solution to the evolu- 
tion equation. 

Proposition 3.3 (Energy estimate). Fix r > 0. Then there exists a constant C 



depending only on f and r, such that for any solution of (32) F we have 



[<7(*)W<Ce-*||<?o|| L 2 + C for t > r. (33) 



Proof. First note that the mild solution of (32) is given by the variation of con- 
stants formula 



g(k,t) = e- H g {k)+ / e- k ^f(-, g(-, s)){k) ds (34) 
for k E {0,1,2,...}. 
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We have, for t > 1, by Plancherel's theorem, 



k=l 

oo 

fc=i 

oo 



< J2k\e- k %(k)\ 2 + Tk / e- fc ( i -)/(-,^(-, S ))(A;)d S 
fc=i k=i Jo 



k=l 

oo 

+£* 

k=l 



" fe( * _s) /(v&(v5))(A:) ds 



f-i 



'C 



fc=i 



e -2fc(t-») d s 
t-i V 



ds 



< Cie-* \\g \\ L 2+C 



\L 2 
oo „ t 



^ /" 

/ /(•>*(•. *))(*) 



ds 



< de 1 « + Q 



(35) 

(36) 

(37) 

(38) 
(39) 

(40) 
(41) 
(42) 

(43) 



If t < 1, the integral from to t — 1 can be disregarded and the integral from 
t — 1 to i runs from to t, with no change in the estimates. In this calculation, 
we have used Holder's inequality in (40) and the fact that j. , e~ fe (* _s ) ds = — | — 
thereafter. □ 

3.2 Existence of pinned and space-time periodic solutions 

First we assert the existence of a stationary solution for zero external driving force. 



Proposition 3.4. Under Hypothesis 3.1, equation (31 )o admits a weak solution 
in H 1 ' 2 . 

Proof. Note first tha if 1 / 2 is compactly embedded in L 2 , independent of dimension. 
One can thus find a minimizer of the energy 

E(g) = \[g?m„ - / f* \{x, 1 )d 1 dx (44) 

z JT n JO 

among functions in He ={wG H 1 / 2 : J Tn u = c} with average c. Now denote by 



G(c) := min E(g) 
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1/2 



(45) 



the energy depending on the fixed average c. Since / has zero average, we have 
G(c + 1) = G(c). Furthermore, G is Lipschitz by a simple comparison argument. 
Therefore, G admits a minimum for some cq € IR. The function g CQ is a weak 
solution to 31), since it minimizes the corresponding energy. □ 



Proposition 3.5. Any weak solution {31) p, F > is classical (and thus also a 
stationary viscosity solution). 

Proof. Plugging the weak solution g of the stationary equation into the variation 
of constants formula for the mild solution of the evolution problem, we find that 
g is a stationary mild solution, since 



-kti 



g( k ) + / e - k ^f{^)){k) ds = e- kt g{k) + / e~ k ^kg(k)ds = g(k) 



for k G {0,1,2,...}. From the regularity properties of mild solutions with initial 
conditions in L 2 , we find that the stationary solution must be classical. □ 

The following theorem asserts the existence of a threshold force, up to which — 
but not above which — a stationary solution exists. 

Theorem 3.6 (Existence of a threshold force). There exists F* > such that 



equation (31) p admits a solution for all F < F* , while it has no solution for 
F > F*. 

Proof. Consider 



$ = {F > such that {31) p has a solution}. 



Clearly, because of Proposition 3.4| $ ^ 0. Also, if F > supv?, then (31)^ has no 
solution. Define, therefore, F* = sup{$} < oo. Two things remain to be shown 
in order to establish the result: 
i) F* e $. 



ii) There is a solution to (31 )p for all F < F*. 

Proof of i): Consider a sequence Fj F* and corresponding solutions gj of (31 ) p. 
such that < f Tn gj < 1. Such a sequence exists, otherwise F* could not be 



the supremum of $. Since {F n } is bounded, Proposition 3.5 asserts compactness 
of {g n } in H 1 . Therefore, there is a converging subsequence whose limit again 
satisfies (|3Tj)f*. 

Proof of ii): Consider < F < F*, and the solutions g* and go to the stationary 
equation with external driving force F* and 0, respectively. By the periodicity of 
(p, and the continuity of the solutions we can assume g* > go. These solutions 



are also super- and subsolutions to (31 )f, respectively. Therefore, there exists a 
solution to (31 1^. □ 
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The next theorem proves the existence of a space-time periodic solution to (32) p 
above the critical force. 



Theorem 3.7 (Existence of a space-time periodic solution). For each F > F* , 



there exists a unique < T(F) < oo and a unique function g(x,t) satisfying (32)p 
such that 

g(x,t + T) = g(x,t) + l. (46) 



Proof. Proposition 3.2 shows global existence of a classical solution. We will split 
the solution up into the evolution of the average and the evolution of the deviation 
of the average and use Schauder's fix point theorem to prove the existence of a 
solution satisfying ( |46| ). 

Define p(t) = j-g(x, t) dx and £(x,t) = g(x,t) — p(t). The functions p and £ 
then satisfy 



P(t) = jf 9 (x,t)dx, (47) 

it(k,t) = -kg(k,t) + fCM^t))(k) fovk>l. (48) 

Consider the initial condition £(x,0) = £o(x) G L 2 , and p(Q) = 0. Since F > F*, 
no stationary solution exists. Therefore, there is T(£ ) > 0, such that p(T(£ )) = 
p(0) + 1, and a constant r, independent of £q, such that T(£ ) > t. This follows 
from the fact that \p\ < \\f\\ Laa . We also have T(£ ) < oo for all £ G L 2 , since 
otherwise one could find a stationary solution to the problem at F > F*, namely 
the pointwise limit for t — > oo of the solution to the evolution problem with initial 
condition £ - 

Now consider the nonlinear operator that advances the solution £ in time, such 
that 

T(£o) = £(-,n£o)). (49) 



and we have 

i(T,k) = e~ kT Uk)+ f T e- k ^f(;g(;sy)(k)d, k > I. (',()) 



From a similar calculation as in the proof of Proposition 3.3 it is clear that 



I|7"(£o)||l2 < e T ||£o|Il2 + C, therefore, for A — C/(l — e T ), this operator maps 



the set ||£o||l2 < A onto itself. The regularity estimate in Proposition 3.3 also 
shows that £(T) is bounded in H 1 / 2 independent of £ > as l° n g as ||Co|lz,a < A. 
The operator T is thus compact and an application of Schauder's fix point the- 
orem yields the existence of a time-space periodic solution which is classical by 



Proposition 3.2 Uniqueness of T and g follow from the Proposition 3.8 below. □ 



15 



Proposition 3.8 (Uniqueness of the space-time periodic solution). The time- 
period T for a solution to (32) p satisfying 

g(x,t + T)=g(x,t) + l (51) 

is unique. Also, the solution itself is unique up to a time-shift so that, given two 
solutions g\ and g<i there exists t such that gi(x,t) = g2(x,t + t ). 

Proof. Assume that there exist two space-time periodic solutions g\ and g2, with 
time constants < T2 < T\ < 00. Since the solutions are continuous and invariant 
under translations by an integer, one can find iV G N such that, for some time 
to, one has «7i(-,£o) < 92(-,to) + N. But since T 2 < T\ there exists a time T 
after which the two solutions would have passed each other, contradicting the 
comparison principle. 

Now, consider a solution G\ with the initial condition <?i(-,to) an d a solution 
G2 with initial condition g^O, ^o)- There have to exist a time T, an integer N, and 
a point Xo such that we have 

Gi(x ,T) = G 2 (xo,t ) + N and Gi(-,T) < G 2 (-,t Q ) + N, (52) 

i.e., the solutions have to touch at some time. Evolving both solutions in time from 
there on, one can see that they touch again after one time period. This, however, 
again contradicts the comparison principle, unless Gi(-,T + 1) = G 2 (-,t). □ 



4 Power laws near the depinning transition 

In this section we investigate the behavior of the average interface velocity near 
the depinning transition. We first note that in the ODE model studied in [TJ, |3], 
namely 

g(t) = <p(g(t)) + F, (53) 

for a Lipschitz continuous, 1-periodic function (p: R — > R and F > F c = — mimp, 
one can obtain any behavior between 

( and v = C\F-F r \\ 



-log|F-F ( 



in leading order for some constant C > 0, where v — (^J^ J is the average 

interface velocity. The limiting case of a logarithmic behavior can be produced by 
taking 

<p{g) = 2 \g- 1/2|-1, 
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Number of Fourier coefficients 1024 

Length over which the interface velocity is averaged 4 

Initial upper bound for F* in bisection 0.5 

Initial lower bound for F* in bisection 

Threshold for accuracy of F* 2 • 10~ 9 

Coefficient of elastic force 0.1 

Time step 1 ■ 10~ 3 

Threshold for stuck interface 1 ■ 10~ u 



Table 1: Parameters used for the numerical examination of the depinning transition 



which yields v = log ^ F -~ ) _l og ^ F _ 1 ^ and thus behaves like _ log \ F _ F | for F close to 
F c = l. 

On the other hand, consider 

for < g < 1/4 
<p(g) = { "I for 1/4 < g < 3/4 

- 3/4) for 3/4 < g < 1, 

which yields v = 1-log( F)+F log( pj+^p^yp log( F _ 1} ■ One can easily see that v = 
2(F — 1) in leading order for F > 1. 

For the non-degenerate case of tp admitting a non-vanishing second derivative 

1 /2 

at its minimum, the square root power-law of v = C \F — F c \ 1 has been shown 




in [U [5] for the ODE case (53). As stated in the introduction, under similar non- 
degeneracy conditions, this power-law behavior has been proved for the case of a 
parabolic model Q by Dirr and Yip in [9]. 

4.1 Numerical method for the investigation of the depin- 
ning transition 

In this section we show some numerical results yielding a good agreement with the 
square-root power law behavior in the case of smooth, non-degenerate obstacles. 
As it was shown in the beginning of this section, this behavior is not generic - in 



fact, also for our model any depinning behavior seen in the ODE model (53) can 
be reproduced by simply picking a heterogeneity <p(x,t) independent of x. 

It follows that the depinning behavior depends very sensitively on the dis- 
cretization of the pinning force, since a piecewise linear discretization of a smooth 
obstacle field for example can easily destroy the original power-law behavior. So, 
starting with an initial configuration g(0) =0 and a fixed applied load F, we nu- 



merically integrate equation (28) using an explicit first-order Euler scheme. This 



scheme is the most appropriate here, since a very small time step has to be chosen 
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Size of pinning sites Critical applied force F* 

1/8 0.0204 

1/16 0.0307 

1/32 0.0150 

1/64 0.0143 

1/128 0.0119 

Table 2: Dependence of the critical applied force on the size of the pinning sites 



in order not to 'jump over' critically pinned states near the depinning transition, 
which immediately renders implicit schemes useless. We also avoid using a higher- 
order scheme, since we want to sample the pinning force at short intervals for the 
numerical integration and not approximate it by a higher-order polynomial. 

The elastic force in our scheme, however, is calculated to high accuracy using 
discrete Fourier transforms. In order to obtain a smooth discretization of the 
pinnign fore, the heterogeneity is constructed using a cubic B-spline. Once the 
interface has traveled a certain length on average (and never got stuck on the 
way), the final time is recorded. This way, a relation between the average velocity 
v and F is obtained. The interface is considered stuck if the L 2 norm of the driving 
force / drops below a certain threshold. This 'inner loop' is repeated with F chosen 
each time through a bisection algorithm, thus giving new upper and lower bounds 
for the critical F* at each run. The program terminates after a certain accuracy 
for determining F* has been reached. In Table [T] the standard parameters for the 
simulation can be found. 



4.2 Simulations 

Experiment 1: General depinning behavior. As a standard example, we 
use a local driving force <p{x\,X2) = ^$(^1,^2), where $ is a potential 
that has smooth dips of a fixed depth and radius at random points. We 
approximate (p by a cubic C 2 spline curve. The exact force used in this 

The constants used for this simulation 



simulation is depicted in Figure 3 (a 



are shown in Table [TJ The evolution of the interface through one period is 
shown in Figure 3(b)| where one can see that the interface spends most of 



its time near the critical stuck state depicted in Figure 3(c) In Figure 3(d) 



the relation between the average velocity, compared to a square-root power 
law is shown. The fit over almost three decades is excellent and at the very 
high end of the applied force one can see that the velocity turns toward a 
linear dependence on the applied force, as expected. 



Experiment 2: Comparison of the depinning behavior for different sizes 
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(a) Force distribution used to ex- 
amine the general depinning be- 
havior. The 'up then down' bumps 
model the x 2 derivative of an at- 
tractive potential well. 



x™ 0.5 




o 1 1 1 1 1 

0.2 0.4 0.6 0.8 1 



X 1 

(c) Interface stuck at the inclusions 
for F = 0.03 



(b) The evolution of the interface 
through one period for F = 0.04. 
Snapshots were taken at equal time 
intervals, so one can see that the in- 
terface spends most of its time near 
the critical pinned state (see Fig- 
ure [3(c) ) 



10" 




1 °io"'° 10" 8 10" 6 10" 4 10" 2 

F-F* 

(d) Power law of the depinning 
transition 



Figure 3: Experiment 1, the general depinning behavior 
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0.2 



0.4 



0.6 



(a) 1/8 



0.2 



0.4 



0.6 



0.8 



(b) 1/32 





(c) 1/64 



(d) 1/128 



Figure 4: The local pinning forces used in Experiment 2. The maximum and 
minimum force are the same for all simulations, the grid on which the spline is 
discretized, however, is of the indicated size. Pinning sites are then distributed 
randomly with constant probability of occurrence. The pinning force for discretiza- 



tion size can be found in Figure 3(a) where also a relation between the color and 
the magnitude of the force is given. 
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F-F* 



Figure 5: The depinning behavior for different sizes of pinning sites. One can see 
that the onset square root power law behavior is pushed towards lower applied 
forces for smaller (and therefore sharper) pinning sites. 

of pinning sites For this experiment, we use five different sizes pinning sites, 
yielding the force distribution as depicted in Figure |4j All other parameters 
are kept as in Experiment 1. The dependence of the critical external force 
F* is shown in Table [2] The depinning behavior is shown in Figure |5j 
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